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Topologically ordered systems are characterized by topological invariants that are often calculated 
from the momentum space integration of a certain function that represents the curvature of the 
many-body state. The curvature function may be Berry curvature, Berry connection, or other 
quantities depending on the system. Akin to stretching a messy string to reveal the number of 
knots it contains, a scaling procedure is proposed for the curvature function in inversion symmetric 
systems, from which the topological phase transition can be identified from the flow of the driving 
energy parameters that control the topology (hopping, chemical potential, etc.) under scaling. At an 
infinitesimal operation, one obtains the renormalization group (RG) equations for the driving energy 
parameters. A length scale defined from the curvature function near the gap-closing momentum is 
suggested to characterize the scale invariance at critical points and fixed points, and displays a 
universal critical behavior in a variety of systems examined. 

PACS numbers: 64.60.ae, 64.60.F-, 73.20.-r, 74.90.-|-n 


I. INTRODUCTION 

In the great number of topological insulators (TI) 
and topological superconductors (TSC) discovered or 
proposedic— , the transition from topologically trivial 
to nontrivial phases can in general be triggered by 
varying all kinds of energy parameters, e.g., chemical 
potentialiSrJ^, hoppingi^, interface couplingi^^— , etc. 
The canonical way of identifying the topological phase 
transition is to first consult the dimension and symmetry 
clas o^^i^^ of the system to see if topologically nontriv¬ 
ial phases are possible, and if so then seek for the value 
of the driving energy parameter at which the spectral 
gap closes at certain momenta, since topological phase 
transition necessarily involves band inversion. Ideally, 
the topological invariant jumps discretely when the sys¬ 
tem undergoes a topological phase transition, therefore it 
seems rather ambiguous to identify any asymptotic crit¬ 
ical behavior from the topological invariant itself. 

Drawing analogy from the Landau order parameter 
paradigm, a profound question is whether there exists 
a universal scaling scheme that can judge the topolog¬ 
ical phase transitions driven by any energy parameter 
in any dimension and symmetry class, and if the scal¬ 
ing procedure is by any means different from Kadanoff’s 
scaling theory^l since the systems under question possess 
no Landau order parameter. Moreover, it would be of 
tremendous usage if the scaling scheme also renormalizes 
the driving energy parameter. The proposal related to 
entanglement entropy first shed a light on this issue^Sr.^. 
On the other hand, since a great number of TIs and TSCs 
are characterized by the topological invariants that are 
often calculated from the momentum space integration of 
a certain integrand function that represents the curvature 
of the many-body state, it is natural to suspect that a 
scaling scheme that renormalizes the integrand function 
can be used to judge topology. Depending on the sys¬ 
tem, the integrand function may be Berry curvature^^— , 
Berry connection^i, or any other quantity whose integra¬ 


tion gives the topological invariant. For the sake of a 
general discussion, hereafter the integrand function is re¬ 
ferred to as the curvature function, and the topological 
invariant calculated from the integration of the curvature 
function is referred to as the winding number (to draw 
analogy to the number of knots in a closed string, and 
to distinguish it from other means of calculating topo¬ 
logical invariants). If a scaling scheme for the curvature 
function exists and gives a RG flow of the driving energy 
parameter that judges the topological phase transition, 
the following criteria are expected: (I) The fixed point 
in the parameter space of the driving energy parame¬ 
ter has a particular configuration of curvature function 
that is invariant under the scaling procedure. (2) Any 
point in the parameter space flows into some fixed point 
while preserving the winding number during the flow. 
(3) Different fixed points correspond to different winding 
numbers. 

In this article, we propose a scaling procedure that sat¬ 
isfies the aforementioned criteria for inversion-symmetric 
systems. The motivation comes from a simple observa¬ 
tion: To know the number of knots that a messy string 
contains, one can either integrate the curvature along the 
string (integrating the curvature function to get winding 
number), or simply stretch the string until the knots be¬ 
come obvious (the proposed scaling procedure). A single 
equation is proposed to implement this knot-tying prin¬ 
ciple in TI and TSC: Let F(k, T) be the curvature func¬ 
tion at momentum k from which the winding number 
in d-dimension C = f d^k F(k, T) is calculated, with T 
the driving energy parameter that controls the topology. 
Given a T in the parameter space, we seek for a new T' 
that satisfies 

F(ko,r') = F(ko + <5k,r) (1) 

where kp is a high symmetry point, and dk is a small de¬ 
viation away from it satisfying F(ko -I- <5k, T) = F(ko — 
dk, r). As we shall see in the examples below, the oper¬ 
ation of Eq. © in inversion-symmetric TI or TSC cor- 


2 


rectly captures the topological phase transition driven by 
any energy parameter. The fixed point is reached by it¬ 
eratively solving r, and by considering a small deviation 
(5k, one obtains the RG equation for T in the form of 
differential equation. The curvature function gradually 
evolves into a fixed point configuration that is invariant 
under this procedure, analogous to a string with its knots 
tight and cannot be stretched anymore. Moreover, for a 
great variety of models, an asymptotic universal critical 
behavior of the curvature function near the gap-closing 
momenta is revealed despite the system in general has no 
Landau order parameter or any short range correlation. 
The critical behavior may be varified by ultracold atoms 
in optical lattices in the case that the curvature function 
is the Berry curvature. 

II. DEVIATION-REDUCTION MECHANISM 

To explain the mechanism behind Eq. o and demon¬ 
strate its analogy to knot-tying, consider a inversion- 
symmetric TI or TSC defined on a (i-dimensional cubic 
lattice. If the system at T and at the fixed point T f have 
the same topology, then the curvature function at T can 
be expanded by 

F(k,r) = F/(k,r/) + F„(k,r) 


This means the deviation at kg changes under this opera¬ 
tion according to the slope of the deviation at kg -I- 5k. If 
the initial value T gives F„(ko,r) > 0 (< 0), i.e., a posi¬ 
tive (negative) deviation at kg, then Fi,(k, T) must curve 
down (curve up) as moving from kg to kg -1- 5k in order to 
conserve the winding number, so < 

0 (> 0) and consequently |E),(ko,r)| is reduced un¬ 
der this scaling procedure. Thus continuously applying 
Eq. ([T|) makes F„(k, E) approaching zero and F(k, T) ap¬ 
proaching E/(k, T f). In other words, the principle behind 
Eq. ([1]) is that the deviation of the curvature function 
from its fixed point configuration is gradually reduced 
under this scaling procedure. As we shall see below, in 
d = 1 systems, this is synonymous to stretching a messy 
string to reduce its messiness until the knots are tight. 
In general. There could be a situation that F„(ko, T) > 0 
and > 0 at a given T so |E„(kg,r)| 

increases at the beginning of the scaling process, but 
|Eu(kg, r)| must eventually be reduced because the scal¬ 
ing process still leads to dk^ Fy (k, T) | < 0 in order 

to conserve winding number. 

Expanding Eq. ([I} in both dT = F' — F and Sk^ = dl 
yields the RG equation. In the leading order, 

dV ia2.F(k,F)|k=ko 

dl ~ 2 drF{ko,T) ' 


= ^/(k, T/) + X! - I n cos ruiki I (,2) 

mi m2 rrid \i=l / 

where Ff {k,T f) is the fixed point curvature function that 
satisfies F/(ko,F/) = E/(kg -|- 5k,F/), i.e., it is invari¬ 
ant under the operation of Eq. o, and Ft,(k, F) is the 
deviation away from this fixed point configuration at F, 
which has Fourier component Am.r = The 

expansion of F„(k, F) is sound because it must not con¬ 
tribute to the integration of winding number 

c = J d^kF{k,r) = j d^kFf{k,Tf) , ( 3 ) 


such that the system at F and at Fy have the same 
topology. Suppose we choose a particular high sym¬ 
metry point such as kg = (0, tt, 0, 0, ...tt), and a small 
displacement away from it along j-th coordinate 5k = 
(0, 0, ..Skj.., 0, 0). Now we apply the operation of Eq. ([T]) 
by using Eq. @. Expanding around kg gives 

E,(kg,F')-E,(kg,F) 

=-im|5fc|^^...^sgn^ko-^m,r , (4) 

mi m2 m^ 


where sgn^^ ko the sign of 11^=1 cos niiki at the chosen 
kg and at the Fourier component m = (mi,m 2 ---md)- 
Note that Eu(kg,r) is the deviation at kg. Writting 
F„(kg,F') — F„(kg,F) = dF„(kg,F) and introducing a 
pseudo scaling parameter dl = dkj , one obtains 


(iE„(kg,F) 

dl 


\dk^Fy{k,T) 


k=ko+<5k 


( 5 ) 


This explains that the proper scaling parameter should 
be defined as quadratic in the displacement 5fc| = dl, 
whereas the dkj = dl in Eq. m is merely introduced to 
demonstrate the deviation-reduction mechanism. Since 
the RG equation involves only one variable F, it is al¬ 
ways possible to map the RG flow into the motion of 
an overdamped particle in a conservative potentia P^'^^ . 
Equation together with the fact that at critical point 
Fc the curvature function diverges at the gap-closing mo¬ 
mentum which is at a high symmetry point kg if the sys¬ 
tem is inversion-symmetrio^l, indicate that the curvature 
function near the gap-closing point kg can be expanded 
by 


E(kg -I- 5k, F) 


F{ko,T) 

l±^2|Jk|2 


( 7 ) 


As shown below, ^ explicitly characterizes the scale in¬ 
variance at the fixed point or critical point. Note that 
although Eq. 0 is similar to the Ornstein-Zernike form 
of correlation function^, ^ does not represent the cor¬ 
relation length of short range fluctuations which is ob¬ 
viously absent in TI and TSC. Instead, it is a length 
scale defined from the curvature function near the gap¬ 
closing point that signals the scale invariance under the 
proposed scaling scheme. Alternatively, one may use its 
inverse k = 1/^ as a momentum scale to characterize the 
scale invariance, which works equally well. 

Two remarks are made before we move to concrete ex¬ 
amples. Firstly, as in any RG procedure, this scaling 
scheme itself does not give a meaning to the fixed point. 
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i.e., it does not tell us whether the fixed point is topo¬ 
logically trivial or nontrivial, which may nevertheless be 
clarified by direct calculation of Eq. m at Ey or at any E 
that flows to the E/. Secondly, E may increase or decrease 
along the RG flow, but the notion of relevant, irrelevant, 
or marginal coupling constants^ does not seem to apply 
to E, since the topological phase transitions are judged 
by the sign change of the RG equation in Eq. which 
is not directly related to how E changes under the scaling 
procedure. 




FIG. 1: (a)The proposed scaling process, Eq. ([T]), applied 
to the topologically nontrivial phase of the SSH model with 
the choice fco = tt. Given an initial St and the corresponding 
Berry connection dk^k{St) (red line), we find a new St' by de¬ 
manding dkffikiSt) at ko + Sk to be equal to that of dk^k{St') 
(orange line) at fco, as indicated by the dash line. This proce¬ 
dure reduces the deviation away from and leads to the fixed 
point configuration dkipk{Stf) (blue line). The ipk expressed 
as a vector field in the complex space for each configuration 
is indicated by colored arrows, (b) The length scale ^ defined 
in Eq. © as a function of St. (c) The RG flow of St. When 
joining the ipk arrows in (a) at a particular St head to tail 
in sequence to form a string, this scaling procedure resem¬ 
bles stretching the string to reveal whether it has a knot, as 
indicated by color arrows. 


III. APPLICATIONS 
A. Su-Schrieffer-Heeger (SSH) model 

To demonstrate the scaling scheme in d = 1, we start 
from the spinless SSEI modeli^ with periodic boundary 
condition (PBC), described by the Elamiltonian 

H = + 6t)c\.^CBi + (t - + h.c. ( 8 ) 

i 

The topology of this model is judged by the winding 
number C of the operator hk{St) = (t + 6t) + {t — 
St)e~''^ = in the complex space when 

one goes through the entire Brillouin zone (BZ). Defin¬ 
ing qk = hk/\hk\ = the winding number is cal¬ 

culated from the curvature function which is the Berry 


connection^ dk^PkiSt) in this case, 

C = j' dk {qk^dkQk) = ^ ^ dk dk^PkiSt) . (9) 
We now search for a new St' by applying Eq. m 


S^k^k(.St )|/c —/cq — dtkPk(.St^\k—ko-\-Sk 

Writing St' — St = dSt, {Sk)"^ = dl, and using dk<Pk{St) = 
arctan (Imh^/Reft.^) one obtains the RG equation at 
the leading order 


dSt 

IT 

dSt 

If 



if fco = 0 , 


if /co = tt , 


( 11 ) 


both reproduce the correct critical points at Stc = 0 and 
fixed points at Stf = ±t, as indicated by the RG flow 
in Fig.[T] Through directly calculating the winding num¬ 
ber, the Stf = t fixed point is topologically trivial and the 
Stf = —t fixed point is nontrivial. The example shown in 
Fig. [1] (a) (red line) clearly demonstrates the deviation- 
reduction mechanism of Eq. where the initial value 
T = St gives Fy{ko,T) > 0 and dkFy{k,T)\k=ko+Sk < 0 
such that Eu(fco, E) is gradually reduced to zero under the 
operation of Eq. O- The choice of different ko means dif¬ 
ferent Fy{ko,T) to start with, hence the speed of converg¬ 
ing to the fixed point configuration is also different, as 
reflected in the two equations in Eq. dl). Note that the 
phase transition occurs when the gap at fcg = tt closes, 
but the scaling scheme works in this model even if one 
chooses ko = 0 that is not the gap-closing point. 

Alternatively, one can introduce the ratio j = (t — 
St) /{t + St) and discuss the topological phase transition 
upon tuning 7 . Using Eq. o with fco = 0 leads to 




= f{l) 


1 ( 

2 U+i; 


( 12 ) 


which well reproduces the two critical points at 7 = ± 1 . 
Thus the proposed scaling scheme is valid whether the 
system is parametrized by St or 7 . 

From Eq. 0 we obtain the ^ deduced from the gap¬ 
closing point ko = TT 





1/2 


(13) 


As shown in Fig.[TJ ^ = 00 at the critical point St = 0, and 
^ = 0 at the topologically nontrivial fixed point St = —t, 
both signature the scale invariance^. We emphasize that 
the scaling scheme here is a procedure akin to knot-tying, 
and scale-invariance means that the curvature function 
converges to a configuration analogous to a string with 
all its knots tight (curvature function at gap-closing point 
stops changing), as indicated in Fig. [T] (c), both have 
different meaning than those in Kadanoff’s scaling the¬ 
ory. The topologically trivial fixed point St = t has 
^ = Xj'Jf,., meaning that while the amplitude of dk'Pki.dt) 
is approaching zero everywhere, its functional form is ap¬ 
proaching the first harmonic cos k. 
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0 (d+1)/(2d+4)B 


FIG. 2: (color online) (a) Berry curvature of the topologically 
nontrivial phase of d = 2 Chern insulators in the continuum, 
which has the form of Eq. (na, at evenly spaced values of 
0.1 < M < 1.0. All the lines evolve to the fixed point con¬ 
figuration labeled by the blue line under the proposed scaling 
procedure, (b) The length scale ^ defined in Eq. 0 diverges 
at Me = 0 and vanishes at M/ = {d+ l)/(2d -|- 4)B. We set 
B = 1 in these plots, (c) The RG flow of M. 


B. Chern insulators in a continuum 

To demonstrate the feasibility of Eq. ([T]) in higher di¬ 
mensions, we consider spinless Chern insulators in the d- 
dimensional continuum^, the low energy sector of which 
has the generic Dirac form 


H{k) =Y,d^{k)^, 


(14) 


where satisfy the Clifford algebra. The curvature func¬ 
tion from which the winding number is calculated, take 
d = 2 and d = 4 as examples, is 

F oc dadxdbdydc , for d = 2 

F oc e°‘^‘^'^‘^dadxdbdydcdzdddyde , for d = 4 (15) 

which has the generic form 

M + Bk^ 


F{k,M) = 


P -b (M - Bk^f 


(16) 


where the prefactor a depends on symmetry and dimen¬ 
sion but is unimportant for our argument. Using Eq. 
0 with the only high symmetry point ko = (0,0) and 
(5k = (SkxySky), writing the energy parameter to be 
renormalized as M' — M = dM while keeping B constant, 
and defining |(5kp = dl, the leading order RG equation is 


dM 

~dr 


d+1 

2d 


1 

M 


d + 2 


B 


(17) 


which has generically a critical point at Ale = 0, and the 
two fixed points Mf = {d+1)/{2d + 4)5 and Mf = —oo 
assuming B > 0. The length scale calculated from Eq. © 




d 1 n\ ^ 

2M2 “ 


vanishes at the fixed points and diverges at the criti¬ 
cal point. The results for d = 2 spinless Chern insula¬ 
tors, which have di = kx, d 2 = fcy, and d^ = M — Bk^ 
are shown in Fig. [21 where the curvature function is the 
Berry curvature^. In contrast to the knot-tying picture 
in d = 1, the scaling procedure in d = 2 is a process to 
stretch the skyrmion texture of d(k) in Eq. (ITT)) with¬ 
out changing the skyrmion number— until the curvature 
function flattens to second order at the gap-closing point 
(blue line in Fig.|2Ka)). 

Two systems of similar kind are Haldane’s d = 2 
graphene model^ and Kane-Mele model^i^. Consider the 
spinless Haldane model whose expansion around K and 
K' points of the reciprocal space of the hexagonal lattice 
is described by the Hamiltonian^ 

3 

/i(ko -b 6k) = —3t2 cos (j) ± 2^1 {6kyax + 6kx(Jy) 

+ (^M + sin Gz (19) 


where upper sign is for kp = K and the lower sign kg = 
K'. Applying Eqs. (fTKI) and ([T]) yields 


dM _ 3 (3ti/2)^ 

dl A M + 3'/3t2 sin (j) 

^ ^ 3(3ti/2)^ 

2 (M + 3\/3t2 sin (j/) 


( 20 ) 


correctly reproducing the critical points at Ale = 
±3-\/3t2 sinej). 
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FIG. 3: (a) The signs of d 3 (ko) at high symmetry points for 

the four phases of d = 2 Chern insulators, (b) The RG flow of 
M by choosing three different high symmetry points ko with 
the same dk = (5fcj;,0). One sees that choosing a particular 
ko captures the transition caused by gap-closing at this point, 
in accordance with the sign change of d 3 (ko) in (a). 


2M2 


1/2 


( 18 ) 
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C. Chern insulators on a cubic lattice 

To discuss the Chern insulators on a d-dimensional cu¬ 
bic lattice, we consider Eq. (imi in d = 2 and d = 4 with 
components^ 


d = 


sin kx, sin fcj,,.... sin kd, M — 2B{d — cos ki) 


i=l 


( 21 ) 


The 2"^ high symmetry points in the first quartet kg = 
(kox, koy ,..., kod) consist of each k^i being either 0 or tt. 
The curvature function expanded around each ko takes 
the form 


E(ko -f dk,M) 


S(dk) 


(- 1 )^- [{M-ABN^) + E{6k)] 

d + 1 1 

{|dk|2 + [(M-4BiV0-S(dk)]2} ' 

( 22 ) 

\i£0 i£7r / 


where X)iG 0 (SieTr) denotes summation over Ski at which 
koi = 0 (fcoi = TJ")) and N.^ is the number of tt’s in {fcoi}- 
Applying Eq. ([T} with dk = SkgS along one particular 
coordinate s yields 


dM 

~dr 


V 2d J M- ABN^ ^ V d y ’ 


(23) 


and hence the generic hxed points Mf = ABNt^ ± (d -|- 
l)/(2d -I- A)B and Mf = =Foo when choosing a particular 
ko. The length scale defined from the gap-closing point 
satisfies 




d -I-1 


T 


{d + 2)B 


(M - ABN^Y M - ABN, 


1/2 


(24) 


the same generic form of curvature function as the Chern 
insulators^ (with the replacement M —>■ chemical poten¬ 
tial and i? — >■ 1/2 X effective mass), and hence practically 
the same critical points and critical behavior. We remark 
that the SSH model^ and Haldane’s graphene model^ 
have been realized by ultracold atoms in optical lattices, 
where the systems can be driven close to the critical 
point. The predicted critical behavior of ^ can be verified 
in the systems where the curvature function is the Berry 
curvature, which may be measurable by detecting the 
anomalous velocity associated with the Berry curvature^ 
or momentum space interferometry technique a'^^d^ . 

IV. CONCLUSIONS 


In summary, we present a scaling procedure to judge 
topological phase transitions driven by any energy pa¬ 
rameter. The procedure is valid for inversion-symmetric 
models in any dimension and symmetry class provided 
the topological invariant is calculated from the integra¬ 
tion of a certain curvature function. Our formalism 
reveals that the concept of scaling in topologically or¬ 
dered systems falls into a completely different realm than 
Kadanoff’s scaling theory in the Landau order parameter 
paradigm. Based on a simple knot-tying picture, the scal¬ 
ing procedure renormalizes the curvature function while 
keeping the winding number intact, through which the 
RG flow of the driving energy parameter T is obtained, 
as well as the RG equation, Eq. under an infinitesi¬ 
mal operation. In essence, the scaling procedure uses the 
divergence of second derivative of the curvature function 
at the gap-closing momentum to find the critical point, 
and uses the flattening of the second derivative to find 
the fixed point, both achieved under a single operation 
of Eq. ([T|). A length scale defined from the Berry cur¬ 
vature near the gap-closing momentum shows an asymp¬ 
totic universal critical behavior that diverges at the crit¬ 
ical point in first power and vanishes at the nontrivial 
fixed point in square root if it is finite 


The top sign in each ± or =F corresponds to the case when 
the component of ko in the scaling direction dk = SkgS 
is fcos = 0, and the bottom sign is when /cqs = tt. Be¬ 
cause 7V,r takes any integer value from 0 to d, generically 
there are d + A critical points located at Me = ABNt^. 
The phase transition at a particular Me takes place when 
the gap closes^Lil at one set of kg’s that have the same 
Ntt, in accordance with applying Eq. o at these ko’s 
to capture the change of curvature function near them, 
as shown for d = 2 in Fig. [3] The chiral p-wave su¬ 
perconductors in the continuum and in the lattice^^ have 


i oc 


|r-r/|i/2|r-r,|-i ifr^^ioo 
ir-Tcl"^ ifr/ = ±oo 


(25) 


for a variety of inversion-symmetric systems examined. 
Applications to a broader class of models, such as those 
driven by interactions or inversion-asymmetric models, 
will be subject to future investigations. 
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